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LONG-TERM  GOALS 

The  overall  goal  of  the  Assessing  the  Effects  of  Submesoscalc  Ocean  Parameterizations 
(AESOP)  initiative  is  to  produce  a  dynamical  framework  to  evaluate  the  impact  of  submesoscalc 
ocean  parameterizations  on  regional  ocean  model  forecasts.  Our  observations  will  improve 
understanding  of  1 )  internal  tides  and  associated  mixing,  2)  submesoscalc  processes  in  the 
surface  boundary  layer,  and  3)  vertical  velocity  (w)  and  mixing  at  an  upwelling  front. 

OBJECTIVES 

Our  objectives  cover  three  themes:  frontal  dynamics,  the  decay  of  submesoscalc  activity  away 
from  the  coast,  and  internal  tide  generation  and  dissipation. 

•  Frontal  dynamics 

•  What  is  the  distribution  of  w  at  the  front? 

•  What  are  the  main  forcing  terms? 

•  What  effect  does  mixing  have? 

•  Submesoscalc  decay  away  from  the  coast 

•  What  is  the  observed  statistical  description  (c.g.,  temperature  gradient  and  vorticity 
probability  distribution  functions)? 


Approved  for  public  release;  distribution  is  unlimited. 


•  How  docs  it  compare  to  a  high-resolution  regional  model  [Capet  et  al.,  2008]? 

•  Internal  tide  generation  and  dissipation 

•  Where  is  the  internal  tide  generated  in  the  Monterey  Bay  area? 

•  How  and  where  is  it  dissipated? 

•  How  well  do  parametcrizations  describe  internal  tide-driven  mixing? 

APPROACH 

Overview  of  observations.  During  the  AESOP  field  program  near  Monterey  Bay,  our  cruise  from 
29  July-28  August  2006  was  divided  into  three  segments:  a  frontal  survey  from  30  July-6  August, 
an  offshore  line  to  examine  submesoscalc  decay  away  from  the  coast  from  6-10  August,  and  an 
internal  tide  survey  from  10-25  August  (Fig.  1).  ScaSoar  was  towed  behind  a  ship  at  4  m  s‘*  and 
made  sawtooth  profiles  to  depths  of  350-400  m.  We  obtained:  1)  hydrography  from  0-400  m  at 
8-m  vertical  and  ~3-km  horizontal  resolution  from  a  SeaSoar  equipped  with  a  CTD,  fluorometer, 
oxygen  sensor,  and  transmissometcr;  2)  conductivity  microstructurc  sampled  at  2048  Hz  from  a 
Transmitting  Microstructurc  System  (TMS)  mounted  below  the  SeaSoar;  and  3)  currents  to 
depths  of  100,  400,  and  600  m  at  vertical  resolutions  of  4,  8,  and  16  m  using  vessel-mounted 
acoustic  Doppler  current  profilers  (ADCP).  Our  measurements  are  of  similar  horizontal 
resolution  to  regional  models  and  can  be  used  to  validate  numerical  model  output  and  their 
parameterization  schemes.  Conductivity  microstructurc  can  be  used  to  determine  the  Cox 
number,  from  which  the  thermal  diffusivity  (Kt)  and  the  dissipation  of  temperature  gradient 
variance  (x)  can  be  calculated.  The  advantage  of  using  the  ScaSoar  to  obtain  microstructurc 
measurements  is  that  a  large  area  can  be  surveyed  rapidly  and  repeatedly. 

Frontal  surveys.  Two  near-synoptic  surveys  were  made  of  a  mesoscalc  upwelling  front  with  at 
least  9  crossings  of  the  front  at  1 1  km  spacing  (Fig.  2).  While  deepening  of  the  mixed  layer  is 
usually  described  well  by  one-dimensional  mixing  models  (e.g.,  KPP  or  Price- Wcllcr-Pinkcl), 
mixed  layer  rcstratification  is  more  complex  and  involves  lateral  secondary  circulations  and 
strong  vertical  velocities.  An  upwelling  front  was  observed  during  restratification  at  small-scale 
by  Lee  and  D’Asaro  (APL)  using  a  neutral ly-buoyant  Lagrangian  mixed  layer  float,  at  kilometer- 
scale  by  Tri Axus,  and  near-synoptically  over  the  mcsoscale  by  SeaSoar. 

Cross-shore  surveys.  Two  surveys  of  --400  km  in  length  were  made  on-  and  offshore  of  Pt.  Sur 
to  examine  the  decay  of  submesoscalc  variability  away  from  the  coast  (Fig.  1 ).  In  particular  we 
will  make  statistical  comparisons  between  our  observations  and  a  high-resolution  idealized 
California  Current  System  model  [Capet  et  al.,  2008].  Probability  distribution  functions  of 
horizontal  gradients  of  temperature  and  salinity  indicate  how  well  models  resolve  the  observed 
fronts.  Probability  distribution  functions  of  observed  vorticity  are  positively  skewed  due  to  large 
cyclonic  vorticities  in  the  mixed  layer  [Rudnick  2001]. 

Internal  tide  survey.  The  internal  tide  survey  repeatedly  surveyed  three  lines  extending  north  of 
Sur  Platform  and  across  Monterey  Canyon  to  a  submarine  fan  (Fig.  1),  all  of  which  arc  believed 
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to  generate  internal  tides  [Lien  and  Gregg,  2001;  Petruneio  et  al,  1998;  Kunze  et  ah,  2002]. 

These  lines  were  repeated  --20  times  eaeh  and  at  different  phases  of  the  tide  (Fig.  3).  Significant 
mesoscale  variability  was  expected  in  this  region,  which  refracts  internal  tides.  Therefore,  three 
adjacent  lines  were  surveyed  near  R/P  FLIP,  SeaSoar  measurements  are  valuable  for  their  broad 
coverage,  which  complements  the  moored  time  series  at  FLIP  and  rapid  expendable  current 
profiles  made  around  FLIP.  Several  AESOP  models,  ASAP’s  HOPS  model,  and  a  modified 
version  of  the  Princeton  Ocean  Model  (POM)  [Carter,  2009]  include  tides.  Model-data 
comparisons  of  internal  tidal  amplitudes,  energy  density,  and  fluxes  are  made  with  POM. 

WORK  COMPLETED 

Vertical  velocity  at  the  front.  Enrie  Pallas  Sanz  (postdoctoral  researcher,  SIO)  has  completed  a 
manuscript  titled  "Frontal  dynamics  and  analysis  of  w  at  a  front  in  the  California  Current,” 
which  will  be  submitted  to  the  Journal  of  Physical  Oceanography  in  November.  Frontal 
ageostrophie  circulations  are  large  and  contribute  to  restratifieation  (Fig.  4).  Objectively-mapped 
observations  arc  fitted  to  dynamics  via  the  quasigeostrophie  Q  equation  and  the  generalized  Q 
equation,  from  which  we  obtain  the  distribution  of  w  and  the  forcing  terms  including  the  surface- 
intensified  mixing. 

Internal  tide  generation  and  dissipation.  Modelled  internal  tidal  amplitudes  and  phase-averaged 
observations  (which  should  be  mostly  due  to  the  internal  tide,  but  include  any  variability  with  a 
time  scale  less  than  the  survey  duration)  agree  qualitatively  (Figs.  5  and  6).  Estimates  of  energy 
flux  and  energy  density  from  SeaSoar  data  in  the  upper  ocean  arc  made  [Martin  and  Rudniek, 
2006;  Cole  et  al.,  2009].  We  will  compare  observed  mixing  at  the  internal  tidal  beam  with  one¬ 
dimensional  turbulence  closure  schemes  (c.g.,  KPP  or  Mellor-Yamada)  and/or  dissipation 
parameterizations  based  on  shear,  stratification,  and/or  strain  [Gregg,  1989;  Kunze  et  al.,  1990; 
MacKinnon  and  Gregg,  2003]. 

Submesoscale  decay  away  from  the  coast.  Work  is  ongoing  on  this  topic,  but  to  date  most  of  our 
efforts  have  focussed  on  the  previous  two  topics. 

RESULTS 

Frontal  dynamics,  w  is  difficult  to  measure  directly,  but  is  important  for  restratifieation  and 
transport  of  biological  material.  For  this  reason,  we  infer  w  from  the  quasigeostrophie  and 
generalized  forms  of  the  Q  equation  (the  latter  with  and  without  parameterization  of  eddy 
viscosity).  Streamers  of  chlorophyll/particulatcs  distinct  from  the  surface  chlorophyll  maximum 
arc  correlated  with  subduction.  w  extrema  of  10  m  day^  are  found  on  the  light  side  of  the  front  in 
survey  1  as  a  combined  effect  of  the  ageostrophie  frontogenesis  and  downfront  winds  [Thomas 
and  Lee,  2005].  Survey  2  is  characterized  by  strong  ageostrophie  deformation  and  subduetion  of 
a  cold  filament.  The  consequences  of  including  vertical  mixing  in  the  generalized  Q,  equation  arc 
a  <50%  enhancement  of  w. 
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magnitude  the  upper  -50  m  and  a  reduction  below.  This  effect  of  vertical  mixing  increases  with 
downfront  winds  and  is  sensitive  to  parameterizations  of  viscosity. 

Internal  tides.  One  or  more  internal  tidal  beams  were  identified  by  calculating  velocity  variance 
along  -60  cross-canyon  sections  derived  from  -20  repeats  of  3  lines  separated  by  2-4  km  (Fig. 

3).  The  beams  appear  to  be  generated  at  the  Sur  Platform,  the  canyon  rims,  and/or  a  submarine 
fan  (Fig.  5).  By  averaging  the  many  repeats  of  these  lines  at  different  phases  of  the  tide,  the 
signal-to-noisc  ratio  was  increased  and  the  velocity  structure  of  internal  tidal  beams  and 
associated  mixing  was  identified  against  the  background  of  an  energetic  mcsoscalc.  Along  these 
beams,  it  appears  Kj  estimated  from  the  microstructurc  data  (Fig.  7)  is  elevated  compared  to 
background  values  away  from  the  tidal  ray  paths  (Fig.  8). 

IMPACT/APPLICATIONS 

When  making  repeat  transects,  our  microstructure  measurements  from  SeaSoar  can  be  used  to 
identify  sites  of  enhanced  mixing  over  a  broad  area  at  -3-km  resolution. 

RELATED  PROJECTS 

Collaboration  with  other  AESOP  Pis  is  ongoing.  Our  results  will  be  useful  to  D’Asaro  and  Lee’s 
frontal  survey.  We  will  be  comparing  our  results  with  Klymak  (UVic)  and  Pinkcl’s  (SIO)  FLIP 
data.  Our  cross-shore  statistics  will  be  compared  with  Capet  ct  al.’s  [2008]  simulations.  We  arc 
Pis  with  MacKinnon  (SIO)  on  the  IWISE  DRl  which  is  investigating  internal  tide  generation  and 
propagation  in  a  sheared  flow.  Carter  (UH),  Gregg  (APL),  and  Johnston  are  organizing  a  session 
at  Ocean  Sciences  2010  titled  "Monterey  Bay  region:  Physical  processes  and  their  impacts.” 

REFERENCES 

Capet,  X.,  J.  C.  McWilliams,  M.  J.  Molcmakcr,  and  A.  F.  Shchepetkin,  Mcsoscalc  to 
submcsoscale  transition  in  the  California  Current  System.  Part  I:  Flow  structure,  eddy  flux,  and 
observational  tests,  J.  Phys.  Oceanogr.,  38,  29-43,  2008 

Carter,  G.  S.,  Barotropic  and  baroclinic  M2  tides  in  the  Monterey  Bay  region,  J.  Phys.  Oceanogr., 
submitted 

Cole,  S.  T.,  D.  L.  Rudnick,  B.  A.  Hodges,  and  J.  P.  Martin,  Observations  of  tidal  internal  wave 
beams  at  Kauai  Channel,  Hawaii,  J.  Phys.  Oceanogr.,  39,  421-436,  2009 

Gregg,  M.,  Scaling  turbulent  dissipation  in  the  thermocline,  J.  Geophys.  Res.,  94,  9686-9698, 
1989 

Martin,  J.  P,  and  D.  L.  Rudnick,  Inferences  and  observations  of  turbulent  dissipation  and  mixing 
in  the  upper  ocean  at  the  Hawaiian  Ridge,  J.  Phys.  Oceanogr.,  37,  476-494,  2007 


4 


Approved  for  public  release;  distribution  is  unlimited. 


Kunze,  E.,  A.  J.  Williams  III,  and  M.  G.  Briscoe,  Observations  of  shear  and  vertical  stability 
from  a  neutrally  buoyant  float,  7.  Geophys.  Res.,  95,  18  127-18  142,  1990. 

Kunze,  E.,  L.  K.  Roscnfcld,  G.  S.  Carter,  and  M.  C.  Gregg,  Internal  waves  in  Monterey 
Submarine  Canyon,  7.  Phys.  Oceanogr,  32,  1890-1913,  2002 

Lien,  R.-C.,  and  M.  C.  Gregg,  Observations  of  turbulence  in  a  tidal  beam  and  across  a  coastal 
ridge,  J.  Geophys.  Res.,  106  (C3),  4575-4591,  2001 

MacKinnon,  J.  A.,  and  M.  C.  Gregg,  Mixing  on  the  late-summcr  New  England  Shelf-  solibores, 
shear,  and  stratification,  J.  Phys.  Oceanogr,  33,  1476-1492,  2003 

Pctruncio,  E.  T.,  L.  K.  Rosenfeld,  and  J.  D.  Paduan,  Observations  of  the  internal  tide  in  Monterey 
Canyon,  y.  Phys.  Oceanogr,  28,  1873-1903,  1998 

Rudnick,  D.  L.,  On  the  skewness  of  vorticity  in  the  upper  ocean,  Geophys.  Res.  Lett.,  28, 
2045-2048,  2001 


Thomas,  L.  N.,  and  C.  M.  Lee,  Intensification  of  ocean  fronts  by  down-front  winds,  7.  Phys. 
Oceanogr,  35,  1086-1 102,  2005 


125  124.5  124  123.5  123  122.5  122  121.5  121  W 


I 


I 


5000 

4000 

3000  £ 
£ 

2000  S' 

a 

1000 

0 


Figure  L  The  n  early-synoptic  survey  of  the  upwelling  front  was  repeated  twice  with  many 
front  crossings  separated  by  11  km  (red).  The  offshore  survey  extended  ^400  km  offshore  to 
document  the  decay  of  submesoscale  activity.  The  tidal  survey  repeated  3  meridional  lines 
about  20  times  to  cover  possible  refraction  of  the  internal  tide  by  topography  or  mesoscale 

flows. 
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Figure  2.  a)  Sea  surface  temperature  (SST)  from  A  VHRR  shows  cold  upwelled  water  near  the 
coast  and  an  offshore  front  Cruise  track  from  survey  1  is  shown  (dots),  b)  Sea  surface 
chlorophyll  a  shows  2  streamers  of  chlorophyll-rich  upwelled  water  advected  along  the  front 
Cruise  track  from  survey  2  is  shoyvn  (dots).  Satellite  data  courtesy  of  NOAA  Coastwatch, 
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Figure  i.  The  tidal  survey  sampled  each  latitude  of  the  3  north-south  lines  at  a  different  point 
in  the  tidal  phase  (SeaSoar  data  in  blue,  TMS  data  in  red).  Since  the  surveys  do  not  resolve 
the  tide,  data  are  averaged  across  tidal  phase. 
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Figure  4.  Vertical  velocity  (w)  from  survey  1  using  the  generalized  Q  equation  at  depths  of  a) 
25.5  m,  b)  56.5  m,  and  c)  192.5  m.  d)  In  a  vertical  slice  along  36.2^N,  upwelling  is  found  on 
the  light  side,  while  downwelling  is  found  underneath  the  sloping  isopycnals  of  the  front.  The 
largest  velocities  are  in  the  upper  150  m  with  velocity  extrema  reaching  10^"^  m  s~L 
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Figure  5,  Velocity  and  displacement  variance  from  each  of  the  about  20  repeats  of  the  3 
north-south  sections  are  scaled  to  provide  an  estimate  of  energy  density.  Topography  is 
shaded  in  grey  at  100  m  increments.  Internal  wave  characteristics  at  the  M2  frequency  are 
shown  (black  lines).  The  zonal  axis  is  elongated.  Based  on  similar  energy  flux  plots,  in  the 
centre  section  one  beam  appears  to  go  upward  and  northward  from  Sur  Platform,  while  in  the 
eastern  section  another  beam  emanates  from  the  north  side  of  the  canyon  and  goes  southward 
and  upward.  The  inset  figure  shows  the  location  of  the  repeated  tidal  survey  tracks  going 

north  from  Sur  Platform. 
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Figure  6.  A  more  complicated  pattern  emerges  from  the  mode!  data  [Carter^  20091,  but  the 
northward  beam  in  the  centre  section  and  the  southward  beam  in  the  eastern  section  are 
confirmed.  Mode!  and  data  show  some  qualitative  agreement,  but  the  magnitude  of  the  energy 
density  is  a  factor  of  3  higher  than  in  the  model. 
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Figure  7.  Microcondiictivity  is  converted  to  temperature  and  gradient  spectra  (blue)  are  fit  to 
the  sum  (red)  of  the  Batchelor  spectrum  (magenta)  and  a  noise  spectrum  (green).  The 
Batchelor  wavenumber  (black  dashed  line)  and  95%  confidence  intervals  (grey  shading)  are 
indicated.  Four  examples  of  these  fits  are  shown. 
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Figure  8.  Median  averaged  thermal  diffnsivity  (Kt)  and  turbulent  dissipation  (e)  are 
calculated  from  approximately  20  north-south  sections  of  the  internal  tidal  survey.  Both  are 
high  along  M2  tidal  ray  paths  (black  lines,  same  as  Figs.  5-6)  emanating  from  the  Sur 
Platform  and  the  northern  edge  of  the  Monterey  Canyon.  Mixing  along  the  ray  paths  appears 
larger,  while  a  low  diffusivities  are  found  away  from  the  ray  paths.  Measurements  in  the 

upper  water  column  are  blanked  out. 
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